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A  Kinetics  Study  of  the  Electron  Cyclotron  Resonance  Plasma 

Oxidation  of  silicon 

J.  Joseph*,  Y.Z.  Hu  and  E.A.  Irene 
Dept,  of  Chemistry,  CB#  3290 
University  of  North  Carolina 
Chapel  Hill,  NC  27599-3290 


Abstract 

The  electron  cyclotron  resonance  plasma  oxidation  of  Silicon 
was  investigated  using  in-situ  during  process  static  spectroscopic 
ellipsometry  and  dynamic  real  time  ellipsometry  at  oxidation 
temperatures  between  80°C  and  400®C  and  at  various  applied  bias'. 
Successful  optical  modeling  of  the  ellipsometric  data  was 
accomplished  using  a  two  layer  model,  in  which  the  top  layer  is  a 
pure  silicon  dioxide  film  over  an  interface  layer.  The  kinetics 
results  are  compatible  with  the  Cabrera-Mott  theory  for  the 
oxidation  by  charged  species  in  the  limit  of  low  electric  field. 
The  effect  of  applied  bias  suggests  that  the  oxidizing  species  is 
O'.  The  energy  activation  is  0.18  eV,  substantially  lower  than  the 
thermal  oxidation  value. 
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I ntr oduc t i on 

Present  silicon  technology  requires  low  temperature 
processing.  For  this  purpose  the  electron  cyclotron  resonance, 
ECR,  plasma  appears  to  be  a  promising  processing  tool,  because  the 
high  ionization  ratio  associated  with  ECR,  yet  the  low  energy  ionic 
species  enables  low  processing  temperatures  with  low  dcunage'.  At 
the  present  time  most  ECR  studies  are  with  film  deposit ion^’^  and 
etching^’*,  and  relatively  fewer  studies  of  Si  oxidation  are 
available®’’’*-®’^®.  In  our  previous  study”  we  investigated  the 
damage  associated  with  ECR  oxidation  of  Si,  and  found  that  the 
damage  is  due  to  the  oxidation  reaction  itself  and  not  typical  of 
ion  bombardment  damage,  and  we  reported  bias  effects  on  the 
oxidation  kinetics.  Since  this  technique  enables  the  production  of 
thin  oxides  at  low  temperatures,  the  present  study  is  aimed  at 
understanding  the  mechanism  for  oxidation,  and  includes  ECR  plasma 
oxidation  of  Si  at  temperatures  between  about  SO^C  and  400®C  using 
in-situ  ellipsometry  for  the  characterization.  In  order  to  obtain 
a  reliable  description  of  the  growth  of  the  oxide  layer,  we 
performed  two  complementary  in-situ  ellipsometric  measurements: 
static  spectrcscopic  and  dynamic  real  time  measurements  at  a  fixed 
wavelength.  The  former  technique  enables  a  precise  determination 
of  the  correct  optical  model  and  the  determination  of  model 
parameters  such  as  damage  layer  thickness,  and  film  and  damage 
layer  composition,  and  the  latter  technique  provides  dynamic 
information  on  film  growth  as  a  function  of  oxidation  time.  We 
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found  that  the  sample  bias  and  temperature  are  the  main  parameters 
for  understanding  the  oxidation  kinetics.  Our  previous  study” 
revealed  that  the  ECR  oxidation  produced  a  complex  layer,  and  that 
multilayer  modeling  is  required. 

Experimental  Procedures 

Figure  1  shows  a  diagram  of  the  apparatus  which  consists  of  a 
high  precision  rotating  analyzer  spectroscopic  ellipsometer  system 
with  an  independent  vacuum  process  chamber  equipped  with  a  homemade 
ECR  plasma  source;  this  apparatus  was  previously  discussed”’”. 
The  plasma  conditions  used  were:  300  W  input  power  at  a  frequency 
of  2.45  GHz  ,  oxygen  pressure  of  5x10"*  torr  .  The  distance  from 
the  microwave  cavity  to  the  silicon  sample  wafer  was  20  cm. 

The  spectroscopic  ellipsometer  measurements  were  taken  at  41 
photon  energies  between  2.5  and  4.5  eV.  This  spectral  region  was 
chosen  because  it  contains  the  main  features  of  the  optical 
spectrum  of  silicon,  namely,  the  two  interband  transition  at  3.4 
and  4.2  eV.  At  lower  energy  the  value  of  A  for  bare  Si  with  a  thin 
oxide  is  close  to  180°  and  therefore  difficult  to  measure  with  a 
rotating  analyzer  ellipsometer”  and,  at  higher  energy  the 
sensitivity  of  the  system  decreases  due  to  insufficient  light 
intensity.  For  all  the  spectra  the  polarizer  was  set  at  25°  which 
optimized  the  accuracy  throughout  the  spectrum.  For  each  sample 
the  offsets  oi  analyzer  and  polarizer  were  redetermined.  This  was 
necessary  to  maintain  the  highest  accuracy,  because  a  small 
misalignment  of  the  sample  stage  occurred  with  sample  changes. 
Each  spectroscopic  measurement  takes  about  10  min  and  during  the 
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measurement  the  plasma  was  stopped.  For  the  oxides  grown  above 
room  temperature  the  spectra  were  obtained  at  room  temperature 
after  cooling  the  sample.  This  was  done  because  of  the 
unavailability  of  reliable  optical  data  for  all  the  components 
present  for  the  higher  temperatures.  At  each  wavelength  a  Fourier 
analysis  of  the  periodic  output  signal  yielded  the  ellipsometric 
parameters  A  and  ^  from  which  the  complex  reflection  coefficient 
is  obtained: 


p  =  tan7exp(iA) 

From  a  literature  database  for  the  known  constituents  of  the  film 
and  substrate  p^  is  calculated  and  compared  with  p^  and  as  a 
figure  of  merit  for  comparison,  an  unbiased  estimator,  5,  is 
calculated  from  the  relationship**: 


6  = 


1 


,1 


2|p  - 

N-P-l  •*** 


where  N  is  the  number  of  wavelengths  sampled,  and  P  the  niimber  of 
unknown  parameters.  A  minimizing  procedure  gives  the  best  fit 
parameters  which  are  film  thicknesses  and  percents  for  the 
constituents  at  the  90  %  confidence  interval.  Throughout  this  work 
the  unbiased  estimators  are  in  the  range  0.003-0.008.  A  value 
outside  this  range  was  considered  as  an  indication  of  problems  in 
the  measurement  or  a  discrepancy  in  the  modelling.  The  real  time 
single  wavelength  measurements  were  carried  out  at  340  nm  (3.65eV)  . 
This  energy  was  chosen  for  several  reasons.  First,  the  values  of 
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A  and  SE^  for  the  bare  silicon  and  thin  oxides  are  in  a  region  where 
the  accuracy  of  the  rotating  analyzer  ellipsometer  is  good. 
Second,  the  dependence  of  the  dielectric  function  of  silicon  on 
sample  temperature  is  the  smallest.  Figure  2  shows  our 
spectroscopic  ellipsometric  measurements  for  a  bare  silicon  wafer 
after  chemical  cleaning  at  three  temperatures,  and  these  data  agree 
with  previous  results  **.  It  is  seen  that  the  spectra  are 
coincident  at  about  3.65  eV.  This  discovery  enables  the  attainment 
of  highly  accurate  ellipsometric  data  despite  temperature  changes 
that  occur  during  processing.  For  sample  temperatures  above 
ambient  the  heating  is  regulated;  but  for  measurement  without 
applied  heating  (refered  to  here  as  room  temperature)  the 
temperature  floats  to  about  80®  C  depending  on  the  duration  of  the 
oxidation  experiments.  With  the  use  of  the  3.65  eV  photon  energy 
we  are  able  to  ignore  the  temperature  changes  in  the  optical 
properties  of  Si  and  obtain  the  changing  thickness  of  the  oxide 
layer.  Third,  there  is  high  sensitivity  to  the  thickness  variation 
of  the  oxide  layer.  Figure  3  shows  a  typical  variation  of  the 
pseudo  dielectric  function  for  different  oxidation  times,  and  at 
3.65  eV  a  proportional  change  is  seen  with  increasing  thickness. 
Lastly,  at  this  energy  the  transparency  of  the  oxide  is  high  and  Si 
is  low  rendering  the  two  major  components  optically  dissimilar 
hence  improving  the  measurement  contrast. 

Oxidations  were  performed  using  about  10  ohm  cm  p-Si  (100) 
oriented  single  crystal  Si  wafers  that  were  cleaned  using  the  RCA 
procedure‘s,  and  immediately  placed  in  the  load  lock  of  the  vacuum 
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chamber.  For  all  samples  an  initial  spectrum  was  obtained  prior  to 
ECR  processing,  in  order  to  insure  the  same  starting  surface.  The 
initial  spectra  have  shown  that  the  cleaning  procedure  is 
reproducible  and  that  the  cleaned  samples  could  be  precisely 
described  by  a  one  film  model  using  the  Bruggeman  effective  medium 
approximation**,  BEMA,  to  be  discussed  below,  and  with  the  film 
composed  of  equal  percentages  of  SiOj  and  amorphous  silicon  and  with 
a  thickness  close  to  1  nm  on  a  crystalline  Si  substrate.  We 
believe  that  the  necessity  to  model  the  native  film  as  a  mixture  of 
oxide  and  a-Si  is  an  indication  of  some  roughness  and/or 
contamination . 
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Optical  Modelling 

For  the  analysis  of  the  ellipsometric  measurements  different 
optical  models  were  used  for  the  spectroscopic  data  and  the  single 
wavelength  data.  This  was  done  for  ease  of  extracting  the  desired 
parameters . 

Spectroscopic  ellipsometry  has  been  shown  to  yield  reliable 
information  for  complex  films  with  the  use  of  stratified  layer 
models*®.  Each  of  the  layers  is  either  a  medium  for  which  the 
optical  properties  are  known,  such  as  SiOj,  or  where  the  properties 
can  be  described  by  a  mixture  of  components.  For  this  latter 
case,  the  optical  properties  are  determined  from  the  dielectric 
functions  of  the  components  and  the  volume  fraction  using  the 
Bruggeman  effective  medium  approximation,  BEMA.  The  BEMA  assumes 
a  discrete  but  compositionally  inhomogeneous  film  where  the  state 
of  aggregation  of  the  components  is  smaller  that  the  wavelength  of 
light  but  not  atomic  points. 

In  order  to  limit  the  number  of  parameters,  we  assumed  that 
the  entire  oxide  o>^er layer  can  be  described  by  no  more  than  two 
discrete  layers  and  that  each  layer  was  composed  of  no  more  than 
two  components.  Therefore,  the  unknown  parameters  in  the  model  was 
limited  to  four;  two  thicknesses  and  two  volume  fractions.  The 
layers  were  assumed  to  be  composed  of  one  or  two  of  following 
components  with  well  known  dielectric  functions:  c-Si“,  a-Si^‘, 
SiOj^  and  voids. 

Four  sensible  models  shown  in  Figure  4  were  evaluated  and  the 
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results  are  summarized  in  Table  I.  From  column  2  which  shows  the 
quality  of  the  fit,  in  terms  of  the  unbiased  estimator,  S,  it 
appears  that  model  4,  made  up  of  a  pure  SiOj  topmost  layer  and  an 
interface  layer  composed  of  SiOj  and  a-Si,  always  yields  the  best 
fit  in  terms  of  the  lowest  unbiased  estimator,  S.  This  table  also 
shows  the  errors  associated  with  the  parameters.  These  values  are 
typical  of  the  results  obtained  throughout  this  work.  Sometimes 
voids  added  to  the  topmost  SiOj  layer  made  a  small  improvement, 
however  always  with  a  volxime  fraction  not  larger  than  the 
associated  error.  For  simplicity  this  fact  was  ignored  and  pure 
SiOj  was  used  for  the  topmost  layer.  In  all  of  our  results  the 
volume  fractions  of  SiOj  and  a-Si  in  the  interface  layer  were  close 
to  50  %.  Thus  the  relevant  parameters  of  this  model  are  the 
thickness  of  the  top  SiOj  film,  L„,,  and  the  thickness  of  the 
interface  layer,  For  some  few  oxidation  parameters  significant 
deviations  to  this  model  were  observed  and  these  will  be  mentioned 
later. 

In  order  to  interpret  the  single  wavelength  measurements  in 
terms  of  oxidation  kinetics  (SiOj  thickness  versus  oxidation  time) 
a  trajectory  method  was  used“.  In  this  procedure  a  measured 
(pseudo)  refractive  index  for  the  starting  cleaned  Si  surface  was 
used  and  the  complex  multilayer  system  is  simulated  by  one  layer 
with  an  adjusted  composite  refractive  index.  For  this  calculation 
the  index  of  the  substrate  is  calculated  from  the  first 
experimental  point,  and  using  this  value  a  trajectory  in  the  A,^ 
plane  corresponding  to  the  growth  of  the  layer  with  an  arbitrary 
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index  is  compared  to  the  experimental  points.  For  this  comparison 
an  error  function  is  defined  as  the  sum  of  the  distance  between 
the  experimental  points  and  the  calculated  trajectory  in  the  A,^ 
plane.  With  the  use  of  a  minimization  procedure  the  index  of  the 
layer  is  varied  to  obtain  the  minimum  of  the  error  function.  The 
thickness  of  the  layer  corresponding  to  each  experimental  point  is 
assumed  to  be  that  of  the  closest  point  to  the  trajectory.  From 
this  procedure  a  calculated  value  for  the  index  of  the  layer  is 
also  obtained.  This  value  could  give  some  information  about 
changes  in  the  layer  during  oxidation,  but  it  is  not  as  precise  as 
the  spectroscopic  results.  We  have  verified  that  there  is  no 
measurable  difference  for  the  growth  of  a  thin  oxide  between  this 
one  layer  model  using  an  effective  index  for  the  substrate,  and  the 
more  exact  two  layer  model,  where  the  substrate  is  described  by  a 
native  oxide  layer  over  bare  c-Si.  It  is  worth  noting  that  the 
thickness  of  the  layer  in  the  one  layer  model  is  only  the  thickness 
above  the  initial  layer,  but  in  the  two  layer  model  the  thicknesses 
of  the  layers  are  measured  from  the  bare  ideal  silicon. 

Experimental  Results 

Oxidation  at  Room  Temperature  for  Different  Sample  Bias'. 

We  found  that  an  applied  (relative  to  ground)  sample  potential 
of  -30V  halted  film  growth  and  the  plasma  became  unstable  above 
+60V,  hence  these  values  provide  the  limits  for  this  investigation. 
The  floating  potential  of  the  ECR  plasma  with  no  applied  bias  puts 
the  sample  at  about  -20V  with  respect  to  the  plasma  potential.  By 
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keeping  oxidation  experiments  under  30  min  the  maximum  sample 
temperature  was  less  than  80"C  for  the  room  temperature  mesurements. 
Figure  5  shows,  the  evolution  of  the  ellipsometric  measurements  at 
different  bias'  in  the  A, ^ plane  during  the  first  20  minutes  of  the 
room  temperature  oxidation.  The  results  for  60V  bias  have  not 
been  plotted,  because  the  trajectory  is  too  close  to  the  30  V 
trajectory.  For  ease  of  comparison,  the  experimental  trajectories 
have  been  shifted  a  small  amount  (the  shift  is  within  the 
experimental  uncertainty  of  the  measurement)  to  coincide  at  the 
same  starting  point.  In  this  Figure  the  solid  lines  correspond  to 
a  calculation  performed  with  a  one  layer  model  using  the  result  of 
the  trajectory  method  for  the  optical  index  of  the  layer.  The 
evolution  of  the  oxide  film  thicknesses  as  a  function  of  the  time 
for  different  bias  values  are  shown  in  Figure  6.  From  these 
results  and  also  from  spectroscopic  measurements  that  follow, 
performed  during  the  very  beginning  of  the  oxidation,  it  appears 
that  two  different  oxidation  regimes  occur  during  plasma 
oxidation. 

The  first  regime  corresponds  to  about  the  first  five  degrees 
shift  of  A  downwards,  or  about  2nm  in  oxide  thickness.  The  key 
feature  of  this  early  oxidation  regime  is  that  the  extent  of  the 
oxidation  is  virtually  bias  independent.  The  bias  does,  however, 
affect  the  duration  of  this  regime:  20s  duration  with  negative  bias 
to  90s  with  a  positive  bias.  The  time  dependence  suggests  that  ion 
bombardment  is  the  main  effect  during  this  period.  With  the  high 
concentration  of  positive  ions  in  the  ECR  plasma,  the  negative 
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sample  bias  would  enhance  the  positive  ion  energy,  and  therefore 
the  importance  of  ion  bombardment.  The  effect  of  the  bombardment 
can  also  be  deduced  from  the  difference  between  the  experimental 
points  and  the  theoretical  calculation  in  Figure  5  where  it  is  seen 
that  this  difference  is  greater  for  the  0  and  the  -30  V  bias. 
Spectroscopic  measurements  for  the  positive  bias  process  at 
intermediate  times  revealed  that  the  best  model  for  the  initial 
regime  is  a  one  layer  model  with  an  equal  fraction  of  a-Si  and  SiOj, 
i.e.  the  same  as  for  the  initial  layer  after  cleaning. 

The  second  oxidation  regime  corresponds  mainly  to  the  growth 
of  the  SiOj  film  as  obtained  from  the  spectroscopic  measurements. 
During  this  regime  the  experimental  points  are  very  close  to  the 
theoretical  trajectory  for  all  bias'.  The  different  solid  line 
trajectories  in  Figure  5  correspond  to  the  use  of  different 
calculated  SiOj  indices.  The  real  part  of  the  refractive  index,  n, 
is  close  to  the  value  usually  given  for  thermal  SiOj,  but  there  is 
a  small  non  zero  imaginary  part,  k.  This  imaginary  part  increases 
as  the  bias  decreases,  and  thus  provides  a  clockwise  rotation  of 
the  experimental  data  in  Figure  5.  This  result  could  be  explained 
either  by  a  concomitant  increase  of  the  thickness  of  the  interface 
layer  and/or  by  some  inclusion  of  unoxidized  species  in  the  oxide 
layer  during  the  growth.  The  spectrosconic  ellipsometry 
measurements  are  concordant  and  show  that  there  is  an  increase  of 
the  concentration  of  unoxidized  optically  absorbing  species  during 
growth  in  the  layer  that  decreases  with  positive  bias.  The  volume 
fraction  of  a-Si  in  the  interface  layer  is  close  to  50%  for 
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positive  or  the  zero  bias,  but  for  the  negative  bias  experiments  it 
was  close  to  75%  in  concordance  with  the  higher  k  value  obtained 
with  the  trajectory  method. 

The  evolution  of  the  thicknesses  as  a  function  of  the 
oxidation  time,  obtained  using  the  trajectory  method,  during  the 
first  20  min  is  shown  in  Figure  6.  The  spectroscopic  measurements 
were  used  to  obtain  the  evolution  of  the  thicknesses  during  the 
second  regime.  For  all  these  spectra  Model  4  always  yielded  the 
best  fit  .  The  results  of  this  analysis  are  plotted  in  Fig  7.  As 
previously  mentionned  the  volume  fraction  of  amorphous  silicon  in 
the  interface  layer  was  found  close  to  50%  for  all  these  spectra. 
The  thicknesses  of  the  interface  layer  are  also  constant  and  nearly 
independant  of  bias.  The  deviations  are  withhin  the  90%  confidence 
limits  which  is  about  0.5  nm. 

Oxidation  at  Different  Temperatures 

With  a  bias  of  +30V,  corresponding  to  the  middle  of  the  useful 
bias  range,  ECR  oxidations  were  performed  at  200‘’C,  SOO^C  and  400°C. 
The  results  from  the  single  wavelength  kinetic  measurements  are 
shown  in  Figure  8.  From  these  results  it  is  clear  that  the  shape 
of  the  growth  rate  is  similar  for  all  the  temperatures;  and  the 
results  of  the  spectroscopic  measurements  are  also  similar  to  the 
room  temperature  results.  However,  there  is  a  difference  in  the 
results  of  the  modeling  of  the  spectroscopic  measurements.  For  the 
longest  oxidation  times  yielding  oxides  thicker  than  25nm  at  300"C 
and  400‘’C,  voids  had  to  be  added  to  the  oxide  layer  to  obtain 
acceptable  unbiased  estimator  values.  This  result  points  out  that 
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even  our  best  model  has  limitations. 

Discussion  of  Results. 

For  thermal  oxidation  of  Si  there  exist  significant  evidence 
that  the  linear  parabolic  oxidation  model  adequately  describes  the 
appropriate  physics  of  oxidation  and  yields  a  reasonable 
approximation  to  the  kinetics^"*.  In  this  model  there  are  two 
oxidation  regimes:  first  in  the  early  regime  the  oxide  growth  is 
limited  by  the  reaction  at  the  Si  surface,  and  later  by  the 
diffusion  of  the  oxidizing  species  through  the  oxide.  In  the  ECR 
plasma  oxidation  at  room  temperature  the  moving  species  are  most 
likely  different  from  those  of  thermal  oxidation  and  the  transport 
of  these  species  is  not  only  due  to  Fickian  diffusion,  but  more 
importantly  to  the  drift  of  charged  species  in  the  electric  field. 
Many  of  these  kind  of  models  also  yield  a  linear  parabolic  rate 
law“  and,  as  a  consistency  check.  Figure  9  shows  a  representative 
plot  of  the  oxide  thickness  vs  the  square  root  of  the  time  for  one 
of  our  thicker  oxide  data  sets.  It  is  seen  that  after  the  first 
few  data  points,  the  results  fit  this  time  dependence.  Thus,  all 
our  results  are  interpreted  with  a  simple  relationship: 


L^=Bt*C 

Where  L  and  t  are  oxide  thickness  and  oxidation  time,  respectively 
and  B  and  C  constants.  Since  all  the  curve  fitting  procedures 
using  this  equation  give  an  insignificant  value  for  C,  only  B  was 
used  for  the  comparison  with  the  models.  Therefore,  using  the 
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total  film  thickness,  Li  =  as  the  thickness  of  the  layer, 
Figures  10  and  11  show  the  kinetics  data  (the  points  are  only  for 
the  static  spectroscopic  data)  and  the  linear  fits  which  provide 
the  parabolic  coefficient  B.  In  Figure  11  the  -30V  bias  curve  does 
not  follow  the  square  root  relationship  well.  Both  the  single 
wavelength  and  the  spectroscopic  measurements  show  that  after  a 
fast  initial  regime,  the  thickness  remains  nearly  constant.  It  was 
also  shown  above  that  for  this  bias  the  layer  is  far  from  a 
transparent  oxide  layer,  hence  the  mechanism  is  not  simply  oxide 
growth  and  the  kinetic  results  are  not  useful  for  comparison  with 
the  other  results. 

According  the  Cabrera-Mott  theory  for  oxidation  by  charged 
species  in  the  limit  of  the  low  field,  the  growth  rate  law^‘  is 
given  by; 


where  L  is  the  thickness  ,  t  the  oxidation  time  ,Ci  and  Cj  are 
constants,  is  the  potential  drop  across  the  oxide  film,  4  is  the 
mobility  which  is  related  to  the  diffusion  coefficient  D  by  the 
Einstein  relationship: 


ZeD=\ikT 

Where  Ze  is  the  electric  charge  of  the  ion,  T  the  temperature  and 
k  the  Boltzmann  constant.  The  diffusion  coefficient  also  depends  on 
the  temperature  T  and  the  activation  energy  E,  by  a  classical 
Arrhenius  relationship: 
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Using  all  these  relationships  the  growth  law  is: 

where  Cj'  is  a  new  constant.  In  order  to  compare  our  experimental 
results  to  this  model,  the  dependence  of  the  growth  rate  with  the 
temperature  and  the  bias  is  examined.  This  theory  predicts  that  the 
product  of  the  temperature  and  the  parabolic  coefficient  follows  an 
Arrhenius  expression,  and  this  is  confirmed  for  our  results  in 
Figure  12.  The  linear  fit  of  these  data  yields  an  activation  energy 
of  0.18  eV.  This  result  which  is  close  to  previous  results  of  0.25 
eV*,  is  notably  lower  that  the  activation  energy  for  thermal 
oxidation  (between  1.4  and  2.3  eV^)  and  clearly  indicates  a  very 
different  oxidation  mechanism.  Figure  13  shows  that  the  parabolic 
coefficient  is  approximatively  linearly  dependant  on  the  bias. 
However,  it  should  be  noted  that  is  smaller  than  the  applied 
sample  bias  V,  because  a  significant  part  of  the  potential  is 
dropped  across  the  plasma  sheath.  This  likely  explains  why  even 
with  a  relatively  high  applied  bias,  the  experimental  results  are 
close  to  the  low  field  approximation.  This  explanation  also 
implies  that  the  relationship  between  V  and  Vox  is  linear  at  least 
in  some  range.  The  increase  of  the  growth  rate  with  the  bias 
reveals  that  the  oxidizing  species  are  negative,  but  contrary  to 


16 


electron  density,  the  density  of  negative  ionic  species  is  low  in 
the  plasma.  In  order  to  explain  our  results,  we  consider  a  multi 
step  kinetic  process.  The  plasma  electrons  first  react  by  a 
favored  reaction  with  the  produce  the  molecular  ion,  ,  that 
is  less  stable  than  O2,  hence  provides  the  atomic  species  0  and  O'. 
O'  can  readily  migrate  through  the  oxide  under  a  positive  sample 
bias. 

From  the  evolution  of  the  interface  layer  during  the  oxidation 
the  best  optical  description  is  as  a  mixture  of  oxide  and  of 
amorphous  silicon.  However,  this  description  is  not  necessarily 
unique  and  one  must  consider  that  several  known  interface 
phenomena  such  as  interface  roughness,  presence  of  suboxides, 
strained  oxide  or  a  combination,  are  all  likely.  From  our 
measurements  we  cannot  decide  which  of  these  are  dominant,  but  from 
the  optical  model  we  know  for  sure  that  this  interface  layer  is 
very  different  from  a  pure  oxide  and  that  its  composition  is  almost 
constant.  The  largest  part  of  this  layer  is  formed  during  the 
early  stage  of  the  ECR  oxidation  and  afterwards  changes  little. 
This  layer  behaves  as  a  marker  in  that  it  neither  moves  toward  the 
surface  nor  spreads  out  in  the  oxide  layer  during  the  oxidation. 
This  strongly  suggests  that  the  ECR  oxidation  occurs  between  this 
layer  and  the  outer  surface  where  the  oxidation  results  from  the 
movement  of  both  the  Si  outward  as  a  cation  in  the  interlayer  and 
the  oxidant  species  inward  as  O'.  In  this  hypothesis  the  interface 
layer  is  the  result  of  the  steady  state  flux  of  Si"^  through  the 
first  layer  of  the  oxide.  That  this  layer  is  not  permeable  to  the 
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oxidant  species  Is  not  readily  explained.  This  model  Is  compatible 
with  the  Cabrera-Mott  theory  which  holds  for  anion  or  cation 
oxidation. 


Conclusions 

1.  We  have  shown  that  the  bias  of  the  sample  Is  an  Important 
parameter  for  ECR  plasma  oxidation.  Positive  bias  Increases  the 
growth  rate  and  decreases  the  Interface  layer  thickness.  The 
damage  resulting  from  the  plasma  occurs  during  the  early  stage  of 
the  oxidation  and  can  be  reduced  with  positive  bias. 

2.  The  oxide  layer  obtained  using  ECR  plasma  exhibits  a  complex 
structure.  From  the  spectroscopic  elllpsometry  measurements  this 
film  Is  composed  of  two  layers:  a  pure  silicon  dioxide  layer  and  an 
Interface  layer.  Optically,  the  Interface  layer  Is  composed  of 
amorphous  silicon  and  oxide  but  futher  work  Is  necessary  to  fully 
elucidate  this  layer. 

3.  The  effect  of  the  bias  and  the  low  activation  energy  (0.18  eV) 
for  ECR  oxidation  definitively  shows  that  the  oxidizing  species  In 
plasma  oxidation  are  different  from  those  of  thermal  oxidation. 
The  atomic  Ion  O'  Is  the  most  probable  oxidant  species. 

4.  The  growth  rate  can  be  modelled  by  a  parabolic  law.  A  parabolic 
relationship  for  the  growth  rate  could  be  obtained  with  different 
models,  but  our  results  show  the  dependence  of  temperature  and  bias 
are  close  to  the  Mott-Cabrera  theory  for  the  oxidation  by  the 
charged  species  In  the  limit  of  the  low  electric  field. 
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List  of  figures 
Figure  1.  Experimental  set-up. 

Figure  2.  Pseudo-dielectric  function  of  Si  at  different 
temperatures.  The  sample  is  a  bare  silicon  wafer  after  chemical 
cleaning. 

Figure  3.  Typical  Si  pseudo-dielectric  functions  at  different  ECR 
oxygen  plasma  exposure  times  with  bias  of  30  V  at  room 
temperature . 

Figure  4.  Optical  models  considered. 

Figure  5.  Evolution  of  the  ellipsometric  A-^  trajectories  during 
the  first  20  minutes  of  oxidation  for  three  bias  and  measured  at  a 
wavelength  of  340  nm.  The  time  between  each  experimental  point  is 
20  seconds.  The  solid  lines  correspond  to  a  model  calculation. 
Figure  6.  Thicknesses  of  the  grown  oxide  layer  vs  ECR  oxidation 
time  for  four  biases  at  room  temperature.  These  curves  are  the 
result  of  data  reduction  at  one  wavelength  using  a  one  layer  model. 
Figure  7.  Results  from  two  layer  modeling  of  in-situ  spectroscopic 
ellipsometry  data  yielding  the  oxide  thickness  and  the  interface 
layer  thickness.  The  measurements  are  for  four  biases  at  room 
temperature.  The  solid  lines  are  only  for  ease  of  viewing.  The 
interface  thicknesses  are  represented  by  ■  ,  a  ,  *  ,  9  , 

respectively,  for  60,  30,  0  and  -30  V  applied  bias'. 

Figure  8.  Thickness  of  the  grown  oxide  layer  vs  ECR  oxidation  time 
for  four  temperatures.  The  bias  is  30  V  and  a  one  layer  model  was 


used. 
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Figure  9.  Thickness  of  the  oxide  layer,  obtained  from  single 
wavelength  ellipsometry  measurements,  as  a  function  of  the  square 
root  of  ECR  oxidation  time,  obtained  at  300°C  for  30  V  applied 
bias. 

Figure  10.  Square  of  oxide  thickness  vs  the  ECR  oxidation  time 
obtained  from  spectroscopic  measurements  for  four  temperatures 
measured  with  30  V  applied  bias.  Total  thickness  is  the  sum  of 
the  oxide  thickness  and  the  interface  thickness.  The  solid  lines 
are  the  result  of  a  linear  fit. 

Figure  11.  Similar  results  as  Figure  10,  but  for  four  applied  bias 
measured  at  room  temperature. 

Figure  12.  Arrhenius  plot  of  product  of  the  temperature  and  the 
parabolic  coefficients  deduced  from  linear  fit  of  the  curves  in 
Figure  10.  The  solid  line  corresponds  to  a  linear  fit  of  the 
experimental  data. 

Figure  13.  Parabolic  coefficients  vs  applied  bias  . 
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Table  I.  Structural  parameters  and  90%  confidence  limits  for  the 
samples  oxidized  with  ECR  oxygen  plasma  at  room  temperature  and 
bias  voltages  of  +30  V  (sample  A)  and  +60  V  (sample  B).  The  models 
are  shown  schematically  in  Fig.4. 
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